Ostrogradski approach for the Regge- Teitelboim type cosmology 
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We present an alternative geometric inspired derivation of the quantum cosmology arising from a 
brane universe in the context of geodetic gravity. We set up the Regge-Teitelboim model to describe 
our universe, and we recover its original dynamics by thinking of such field theory as a second-order 
derivative theory. We refer to an Ostrogradski Hamiltonian formalism to prepare the system to its 
quantization. Our analysis highlights the second-order derivative nature of the RT model and the 
inherited geometrical aspect of the theory. A canonical transformation brings us to the internal 
physical geometry of the theory and induces its quantization straightforwardly. By using the Dirac 
canonical quantization method our approach comprises the management of both first- and second- 
class constraints where the counting of degrees of freedom follows accordingly. At the quantum 
level our Wheeler-De Witt equation agrees with previous results recently found. On these lines, we 
also comment upon the compatibility of our approach with the Hamiltonian approach proposed by 
Davidson and coworkers. 
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I. INTRODUCTION 

The concept of a relativistic extended object as a sur- 
face immersed in a bulk has increased the interest in 
physics due to its wide range of applications. One can 
model, for example, the smallest physical entities, like 
quarks, as vibrations of strings up to the entire universe 
as a relativistic extended object. Along a related line, 
with the advent of brane world universes, cosmology in 
the presence of extra dimensions has been the subject of 
intense research. In fact, the idea that our universe could 
be a 3 + 1 dimensional surface embedded in a higher di- 
mensional spacetime was set up by Regge and Teiltel- 
boim (RT) a long time ago [ij and pursued by many 
authors @, [H, 0, E| . The scope of such a model is that 
gravitation can be described in a point- or stringlike fash- 
ion, as the worldvolume swept out by the motion of a 
three-dimensional spacelike brane evolving in a higher- 
dimensional bulk spacetime Recently, the RT brane 
model has been considered as one of the two main pil- 
lars of a unified branelike theory Q , where the Randall- 
Sundrum brane theory [tJ is included. When one ad- 
dresses this issue in a Minkowski spacetime, the model is 
named geodetic gravity, and it has been extensively stud- 
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ied by Davidson and coworkers [1, 0, [T(| E| ■ Although 
the RT model is not the most popular theory for brane 
world universes (at the end of last century, there was a 
revival of the idea that our universe could be a hypersur- 
face; see, for example, @,[Hj]), it is very stimulating while 
thinking in the spirit of brane gravity a la string. The 
cosmology that arises from this model is interesting in its 
own right since it provides an alternative route to better 
understand classical cosmology in extra dimensions, and 
also it supplies a compelling model to apply the canonical 
quantization methods. Indeed, in the context of quan- 
tum brane cosmology [l3| our universe can be explained 
through a tunneling process where the well-known prob- 
lem of boundary conditions of four-dimensional cosmol- 
ogy is solved [IJ, [l5| . 

In most field theories the action depends usually on 
the fields and their first derivatives. By contrast, the RT 
model is a genuine second-order derivative model in the 
field variables, which are the embedding functions rather 
than the induced metric. Generally, one identifies and 
neglects a surface term associated with the linear depen- 
dence of the accelerations. Similarly as in general relativ- 
ity, it is a well-known fact that a "harmless" surface term 
can be neglected or removed at the beginning as occurs 
with the well-known Gibbons-Hawking- York term into 
the action. Whichever field context, the extremization 
of the corresponding action yields equations of motion of 
second order in derivatives in the field variables. Thus 
the RT model, like the Einstein-Hilbert action for general 
relativity, becomes transformed in an effective first-order 
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field theory. However, by what formally appears to be a 
customary procedure, to follow such a strategy raises im- 
portant limitations especially in the Hamiltonian frame- 
work for the RT model where it leads to certain trou- 
bles, as was noticed first by Regge and Teitelboim, due 
mainly to the fact that the scalar constraint is not writ- 
ten down in a closed form straightforwardly. In pursu- 
ing this endeavor, Davidson and coworkers tackled the 
problem successfully. They considered an extra nondy- 
namical canonical field A in the first-order Hamiltonian 
framework in order to get quadratic constraints of the 
phase space that recuperate the dynamics accordingly 
@) OH EH- The explicit handling of the quantum RT 
model is made possible by extending the ordinary phase 
space, which in turn provides a wealth of information of 
the cosmology that this model possesses. 

In the present paper we consider an alternative for- 
mulation for geodetic gravity which is strongly based in 
the Ostrogradski program for higher-order derivative the- 
ories |16l |17| . For second-order theories this approach 
treats the velocities as independent fields. This is an un- 
conventional viewpoint for the RT theory, and one might 
therefore wonder if such a description is viable at all since 
this does not necessarilly represent a shortcoming: for 
this special case the addition of more degrees of freedom 
is physically more accurate, but it means then that the 
first-order theory is incomplete in some sense. For this 
reason, it seems promising to start directly from the full 
RT model instead of omitting the surface term a priori. 
Hence, we pay close attention to a Hamiltonian approach 
for geodetic gravity constructed by appliying the Ostro- 
gradski scheme which in turn leads to the correct dy- 
namics. In particular, it is of a great interest to use the 
full model straightforwardly for obtaining the quantum 
approach for brane cosmology. Our intention is to cope 
directly with the inherent second-order derivative nature 
of the RT model. As discussed below, we gain certain im- 
provements of clarity by the use of this formalism in com- 
parison with previous works. Contrary to the standard 
quadratic form of the constraints for ordinary first-order 
reparametrization invariant theories, in the Ostrogradski 
approach for RT field theory the constraints are projec- 
tions of the momenta along the velocities as well as along 
the unit spacelike normal vector to the brane. To illus- 
trate our development we specialize our considerations 
to a minisuperspace model where the inherent gauge in- 
variance under the reparametrization of time is evident. 
We show that the canonical Dirac constraint quantiza- 
tion of this model casts into a satisfactory Wheeler-De 
Witt (WDW) equation on the wave function for a brane- 
like universe. The handling of the quantum approach is 
made possible by a canonical transformation which re- 
sults to be a Lorentz rotation in phase space. Such a 
transformation brings our constraints into a physically 
meaningful set which enable us to follow the standard 
Dirac constraint quantization programme. Our quantum 
treatment hence leads to a well-defined Wheeler-DeWitt 
equation which, even though it is technically complicated 



to solve, presents the right behavior for the quantum po- 
tential, estimating the accuracy of our approach. 

The outline of the paper is as follows. In Sec. II we 
briefly review some geometrical aspects of the RT model 
for a general <i-dimensional brane, which are of inter- 
est for the rest of the paper. This section will serve to 
explain our notation and to gain insight into our Hamil- 
tonian approach. In Sec. Ill we adapt our approach 
to a minisuperspace model in which we specialize to the 
geometry generated by the Friedman-Robertson- Walker 
(FRW) model. We explicitly give the Lagrangian density 
associated to the RT model which includes the surface 
term. Sec. IV deals with the Ostrogradski approach for 
the model we are considering, and we develop the corre- 
sponding constraint analysis. In Sec. V we propose the 
gauge-fixing for the model in order to completely identify 
the structure of the reduced phase space. In Sec. VI, we 
study the quantization of our model within the scheme 
of Dirac quantization. Finally, in Sec. VII we draw some 
conclusions. As a general feature, our presentation avoids 
cumbersome notation and is intended to be index-free as 
possible. 



II. REGGE— TEITELBOIM MODEL 

Consider a brane £ of dimension d, evolving in a fixed 
Minkowski N dimensional background spacetime with 
metric rj^. Its trajectory, or worldvolume m of dimen- 
sion d+1, is described by the embedding x M = X^(^ a ), 
where x M are local coordinates for the background space- 
time, £ a local coordinates for m, and A" M the embedding 
functions (fi, v = 0, 1, . . . , N — 1; a,b — 0,1, ... ,d). We 
denote by e M Q = d a X^ the tangent vectors to m. In this 
framework we introduce N — d — 1 unit normal vectors to 
the worldvolume, denoted by n^i (i = 1,2,..., N— d— 1). 
These are defined implicitly by n l ■ e a = 0, and we choose 
to normalize them as • rij = &y . 

The RT model for a d-dimensional brane S is defined 
by the action functional 

Sbt[X} = ^[ d d +^^—gK-f d d +^J—gk, (1) 

where the constant a has dimensions [L\^ d \ g denotes 
the determinant of the induced metric g a b = r/^ e^ae^b — 
e a ■ eb- We have also included in this action a cosmolog- 
ical constant term, A. The extrinsic curvature of m is 
K ab l = -n l ■ D a e b , where D a = e^aD^ and is the 
covariant derivative in the bulk spacetime. The mean 
extrinsic curvature is given by the trace K l — g ab K a b l 
where g ab denotes the inverse of g a b- The scalar cur- 
vature 1Z of m can be obtained either directly from 
the induced metric g a b, or, in terms of the extrinsic 
curvature, via the contracted Gauss-Codazzi equation, 

n = k i k, ~ K^Kf mm 

The response of the action ((T|) to a deformation of the 
surface X — > X + 5X is characterized by a conserved 
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stress tensor [20j 

f a ^ = -( a g ab + Ag ab )e^, (2) 

where Q ab = lZ a b ~ \H-9ab is the worldvolume Einstein 
tensor, with lZ a b being the Ricci tensor. This quantity 
will provide relevant physical information especially with 
consistent conservation laws. Following the line of rea- 
soning of [2(3], the classical brane trajectories can be ob- 
tained from the covariant conservation law, V a / a/J = 0, 
where V a is the covariant derivative compatible with the 
induced metric g a b [H|. This yields [13] 

T ab K ab l = , (3) 

where T ab = uQ ab + Ag ab . In fact, T ab corresponds 
to the intrinsic stress tensor defined in the usual way 
by —2/^/^g(SSnT/Sg a b)- Its conservation is supported 
by the Bianchy identity. The equations of motion ([3]) 
are of second order in derivatives of the embedding func- 
tions because of the presence of the extrinsic curvature. 
This is so even though in the scalar curvature TZ we have 
the presence of the extrinsic curvature. Owing to the 
reparametrization invariance of the RT model, there are 
only D — d — 1 independent equations, along the normals; 
the remaining d + 1 tangential components are satisfied 
identically, as a consequence of the reparametrization in- 
variance of the action ((T|). 

An important quantity constructed with the conserved 
stress tensor is given by 

tt" = V af a M = - { a g ab + Ag ab ) Va e» b (4) 

where rf 1 stands for the timelike unit normal vector to the 
brane E when it is viewed into m [23j. In fact, Eq. (U) is 
nothing but the conserved linear momentum associated 
with the Noether charge of the action ([1]) specialized to 
background translations 20]. The E basis, {e M A, ?7 M , n^i} 
satisfies the completeness relation 

if v = l n\ - rfrf + h AB e» A e v B , (5) 

where Hab — gab£ a A£ b B is the spatial metric on E and 
e a a are the tangent vectors to E, (A, B = 1, 2, . . . , d). 
The vector r?^ stands for a timelike unit vector to E (see 
Refs. [2l|, [24| for more details). 

In presence of other possible matter sources with stress 
tensor = (— 2/ y r ^g)5S m /5g a b, where S m is a mat- 
ter action, we do not expect considerable modifications 
in our approach. The equations of motion ((3|) remain 
unchanged in form. It is sufficient to add the matter 
stress tensor to the original one described in ([3]). Sim- 
ilarly, the conserved linear momentum (j4|) is unaffected 
in form when another type of matter is included. It gets 
an additional contribution of the form 7r^ — —T^rjae^b- 
These nice features allow us to develop straightforwardly 
a Hamiltonian analysis without substantial changes un- 
der the inclusion of matter fields. This fact was also 
noticed in [Til ] . 



III. MINISUPERSPACE MODEL 

We turn now to restrict the RT model itself (Q} to 
the case of a minisuperspace model. Consider a 3-brane 
E, evolving in a 5-dimensional Minkowski spacetime, 
ds 2 = —dt 2 + da 2 + a 2 dfl 2 : where dQ 2 stands for the 
metric of a unit 3-sphere, i.e., dQ\ = dx 2 + sin 2 \d0 2 + 
sin 2 x sin 2 9d(f> 2 . For the sake of simplicity, we choose the 
function sin 2 x i n dQ^ to consider a closed universe. If 

x» =X»{C) = {t{T),a{T), X ,e i 4>) (6) 

is a parametric representation of the trajectory of E, we 
assure that the geometry of the worldvolume generated 
is that of the FRW case. According to the cosmology 
jargon, a(r) is known as the scale factor. 

The basis adapted to the worldvolume is given by the 
four tangent vectors e^ a (a = 0,1,2,3) together with 
the unit spacelike normal vector 

% = i(-<M, 0,0,0) , (7) 

where the dot stands for derivation with respect to r. For 
short in the notation we have introduced the quantity, 
N = y/i 2 — a 2 , which coincides with the lapse function 
when we perform an ADM decomposition of the action 

© pin. 

The metric induced on the worldvolume is given by 

ds 2 = g a bd£ a d£ b = -N 2 dr 2 + a 2 dQ 2 . (8) 

The spatial components of this metric correspond to the 
metric associated to E when this is described by its em- 
bedding in the worldvolume itself. Furthermore, for this 
latter parametrization, we have rj a — l/N(l, 0, 0, 0) such 
that 5a6 e<W = 0. 

The Ricci scalar associated with the metric (]8]) reads 

Tl = -^—r (adi - adi + N 2 i) . (9) 

The linear dependence that the Ricci scalar possesses in 
the accelerations of the variables t(r) and a(r) is partic- 
ularly remarkable. 

The Lagrangian density C = ^—g(jlZ — A) thus be- 
comes 

Tata , ■ ■■ ■■, 9 -x TNa 3 . , . 
C = N3 (aat - adt + t 6 - d 2 t) — A , (10) 

where T = 3sin#sin 2 x. Thus, the RT action special- 
ized to spherical configurations, in terms of an arbitrary 
parameter r, is reduced to 

Srt = 67r 2 J dr L(a,a,'d,i,i) , (11) 

where the Lagrangian function is given by [34| 

L(a, a, d, i, i) = (adi - adi + N 2 i) - Na 3 H 2 , (12) 
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where we have introduced the constant quantity H 2 := 
A/3a. Thus, we have only a(r) and t(r) as independent 
dynamical variables. Despite the acceleration depen- 
dence in the Lagrangian, as characterizes second-order 
derivative theories, the equations of motion remain sec- 
ond order in the field variables (see Eq. |3])). 

We proceed now to evaluate both the Einstein and 
the extrinsic curvature tensors of the worldvolume as de- 
scribed by the metric ([8]). We have the nonvanishing 
explicit components 



such that the highest momentum spacetime vector is 



g\ = 



m 2 



a 2 N 2 

i 2 

a 2 N 4 







2ai i (?) 


■f TV 2 







and 



K 4> 



4> 



N^dr 
i 

aN' 



Clearly, we can read off immediately the spatial compo- 
nents of the extrinsic curvature tensor as well as its mean 
extrinsic curvature given by JC = h AB Kab = 3t/aN. 

As dictated by Eq. there is only one equation of 
motion given by 



d_ 

d^ 



n 2 e 



where we have introduced the functions 



9 := t 2 - 3N 2 a 2 H 2 . 
$ := it 2 - N 2 a 2 H 2 . 



(13) 



(14a) 
(14b) 



for simplicity. Eq. (|13p obviously only involves second 
derivatives of the field variables a and t. For any solution 
for a(r) we have a gauge freedom to choose a function 
for t(r) as we will see below. 



IV. OSTROGRADSKI— HAMILTONIAN 
APPROACH 

A deeper insight of the phase space structure of the 
theory defined by the Lagrangian (fT2"j) is achieved by an 
Ostrogradski procedure for higher-order derivative sys- 
tems. (A complete description of the Hamiltonian formu- 
lation for branes whose action depends on the extrinsic 
curvature of their worldvolume is provided in [2ll . 
The highest conjugate momenta to the velocities {i, d} 
are, respectively, 



Pi 

Pn 



dL 
lit 



a at 

' N 3 



dL a 2 i 2 



dd 



N 3 



(15) 
(16) 



P 



at, ■ % at 

= —(-^,0,0,0) = ^ 



N 2 n ^ 



(17) 



Though this momentum has not a direct mechanical 
meaning it will become important to achieve a Legendre 
transformation in order to obtain the Hamiltonian func- 
tion for our system (see Eq. (|24|) below). Note that the 
momentum is directed normal to the worldvolume. 

The conjugate momenta to the position variables {t, a} 
are, respectively 



Pi 



Pa 



9L _d_ (dL 
di ~ dr\ di 

^[d 2 + N 2 (l-a 2 H 2 )] 

dL d (dL 
da dr \ dd 

-^[d 2 + N 2 (l-a 2 H 2 )] = 



ft, 



(18) 



(19) 



Important to note is the fact that both momenta, pt and 
p a , are from a totally different nature. Indeed, while 
the momentum p t is not influencied at all by the surface 
terms (as expected), the momentum p a is obtained by 
two contributions: one coming from the ordinary theory 
and the other by a surface term. In this way, we can 
denote the momentum p a as 



where 



Pa 



Pa 



aa 
~N3 
2aa 
~N~ ' 



Pa + Pa ■ 



N 2 3 - a 2 H 2 



(20) 

(21) 
(22) 



It is crucial to recognize them as the canonical momen- 
tum worked out in [9j and as the momentum conjugated 
to the a(r)-variable when considering as the Lagrangian 
only the surface term L s = d/dr(a 2 a/N), respectively. 

To see the geometrical structure of this momentum, it 
is convenient to write p^ as 

P» = yH,d, 0,0,0) = jj,. (23) 



We realize that this momentum is identical with the vec- 
tor a ^ defined on Eq. ((4]), which is the projection of the 
conserved stress tensor along the unit timelike normal 
vector rf l to S. 

The appropiate phase space of the system, F := 
{t,a,i,d;pt,Pa, Pt, Pa} , has been identified explicitly. 
Thus in T, the Ostrogradski formalism yields the canon- 
ical Hamiltonian 



p- X + P ■ X - L=p a a+ Pt i + J n , (24) 
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where we have defined 

Jn = — [a 2 + N 2 (1 - a 2 ff 2 )] = ^Cn . (25) 

This potential- like term results an implicit function of the 
phase space variables in the combination N 3 P 2 . At first 
glance, this may look as an unnecessary complication to 
write the phase space quantities in terms of f2, but this 
quantity results in a physical observable: It is nothing but 
the conserved bulk energy. Indeed, squaring the energy 
equation (1181) . results in the so-called evolution master 
equation [H, 0, 0, Q 



N 2 + a z = jN 2 a 2 H 



(26) 



where 7 = 7(a) satisfies the cubic equation 7(7 — l) 2 = 
n 2 /a s H 6 . 



A. Constraint analysis 

Since we are dealing just with a second-order deriva- 
tive theory linear in the accelerations, we have two pri- 
mary constraints given by the definition of the momenta 



itself, and hence 



A - 



= 0. Instead of 



these constraints, here we will follow a different but con- 
venient route. We choose to project the momentum (|17[) 
along the velocity vector as well as the normal vector to 
S, where in general = n^(X v ). This is supported 
by using the geometrical completeness relation ([3]) in 
<t>n = r l^4> v ■ Thus, we get the primary constraints 



Ci = P ■ X = 
C 2 = NP-n- 



a 2 i 



0. 



(27) 
(28) 



Therefore, the total Hamiltonian which generates time 
evolution of the fields is 



Ht — Plq + Ai C\ + A2 C2 



(29) 



where Ai and A2 are Lagrange multipliers enforcing the 
primary constraints. 

As customary, time-evolution for any canonical vari- 
able z £ T reads 



{z,H T } 



(30) 



on the constraint surface, where the generalized Poisson 
bracket for any two functions F(z) and G(z) in F is ap- 
propriately defined as 



{F,G} = 



dF dG dF dG dF dG 
dt p t da p a di P t 



d£dG 
~da~~Pa~ 



(31) 



Important to mention is the fact that the total Hamil- 
tonian leads us directly to the right equations of 



motion ([13]) through the conventional Ostrogradski ap- 
proach for higher-order derivative systems [17] . 

We also note that under the symplectic structure (|31jl , 
the constraints (f2"T)) and result to be in involution, 
{Ci,C2} — 0. According to the Dirac program for con- 
strained systems, both C\ and C 2 must be preserved by 
the evolution which demands the existence of the sec- 
ondary constraints 



C 3 = H = p ■ X + N [ a 3 H 2 - \n 2 P 



p ■ n. 



(32) 
(33) 



The vanishing of the canonical Hamiltonian is expected 
courtesy of the reparametrization invariance of the RT 
model. Hence the canonical Hamiltonian Hq generates 
diffeomorphisms normal to the worldvolumc. The sec- 
ondary constraint (|33[) is characteristic for every brane 
model linear in accelerations. The process of generation 
of further constraints is stopped at this stage since C3 
is preserved under evolution and the requirement of sta- 
tionarity of C4 only determines a restriction on one of the 
Lagrange multipliers, namely, A2 = N 3 £l/a 2 &. Thus, we 
are dealing with a wholly constrained theory with first- 
and second-class constraints, which is a consequence of 
the rich gauge symmetry of the RT model. The distin- 
tive feature of the constraints (|27[) and (|28p instead of 
<j)^ is that C\ and C2 are constraints that naturally gener- 
ate the relationships (|3"!?|) and as befit a higher-order 
derivative brane theory [24| . 

Following Dirac's program, the set of constraints 
should be separated into subsets of first- and second-class 
constraints [26] . It is quite well known that for each pair 
of second-class constraints there is one degree of freedom 
which is not physically important and has to be removed 
from the theory, and for each first-class constraint one 
degree of freedom is removed. For our system we have 
two first-class phase space constraints 



T 2 



Ci, 

N 3 n 

a 2 $ 



■C2 + c 3 



(34) 
(35) 



and two second-class constraints. The selection of the 
second-class constraints is straightforward (see the Ap- 
pendix). We choose them as 



S 2 



c 2 



(36) 
(37) 



Note that as we have two linear independent first- 
class constraints, we have the presence of two 
gauge transformations in the RT model. In the 
Appendix we discuss more thoroughly the Pois- 
son brackets among the phase space constraints. 
The counting of degrees of freedom is as follows: 
dof = [(Total number of canonical variables) — 2 x 
(first — class constraints) — (number of second — class 
constraints)] /2 = [8— (2x2) — 2]/2 = 1, which agrees with 
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the number of normals to the worldvolume [24|. Such a 
single degree of freedom can be identified as the scale 
factor a(r). 



V. GAUGE-FIXING 

According to the conventional Dirac scheme, in or- 
der to extract the physical meaningful phase space for 
a constrained system we need a gauge-fixing prescrip- 
tion which entails the introduction of extra constraints, 
avoiding in this way the gauge freedom generated by con- 
straints (|34|) and ([35| . To achieve this we will consider 
the conventional cosmic gauge condition 



fx 



N-1 = V¥ 



lwO, 



and the generalized evolution Eq. ([26 



2„2 



tp 2 = + af --yN^H'a 



0. 



(38) 



(39) 



where the ss symbol stands for weak equality in the Dirac 
approach for constrained systems [26] . From the geomet- 
ric point of view, this set of gauge conditions is good 
enough since the matrix ({J- ', ^1,2}) is nondegenerate in 
the constraint surface. Indeed, under the Poisson bracket 
structure (|3ip. it is straightforward to show that gauges 
tpi and tp 2 form a second-class algebra with the con- 
straints T\ and T 2 



Wx,?x} = fx +1, 
{fx, ^2} = 0, 

{f2,F 2 } = F(a,a,i) , 



(40) 



where F(a,a,i) is a complicated function [35[. Conse- 
quently, velocities t and a must be discarded as dynami- 
cal degrees of freedom. 

The use of the completeness relation ([5]) results efficient 
at this level: It allows us to express the quantity P 2 
by the equivalent expression — (P • r/) 2 + (P ■ n) 2 . This 
suggests the implementation of the following canonical 
transformation to a new set of phase space variables: 



N := y/i 2 - a 2 



n N -.= -(p.x). 



v := arctanh 



Pa 
Pt 



(41) 



n„ 



N(P ■ n) , 



together with the coordinate transformation Z M := A M — 
{X^,v} Tl v , while the momenta p^ remains unaltered. 
Such transformation can be physically interpreted as a 
Lorentz rotation in phase space which, straightforwardly, 
preserves the structure of the canonical Poisson brackets 



{n,u n } = 1 = {v,n v } , 

{Z», Pu } = 5\, 



(42) 



as expected. A transformation of the same kind was con- 
sidered in [25| for a quantum treatment of a kink model. 
To understand in a more transparent way the geomet- 
ric content of this transformation let us go through the 
intrinsic angular momentum density of the RT model 
calculation carefully. In terms of the new phase space 
variables the intrinsic angular momentum is given by 
= VaM^ v = NP^X^ 0. Such momentum is 
conserved in the sense that V a A4'* flL ' — 0. This con- 
tribution shows up as an effect of the finite width of 
the brane in comparison with flimsy branes described by 
Dirac-Nambu-Goto action. The only nonvanishing con- 
tribution to the intrinsic angular momentum reads 



Mi m = 7; NU v , 



(43) 



which in turn tell us that our new variable II„ can be 
thought of as an angular momentum component. 

In attempting to take into account the new phase space 
variables, the first- and second-class constraints ()34II37|) 
become 



T x = An 



N 



(44) 



;a = i> x - .v ( a 3 H 2 + -^N 2 n 2 N - -^ n 2 , ) . 1 4r, , 



and 



Si 



n, 



N 



0, 



(46) 
(47) 



... 2aa ■ 

o 2 = p a t + pta = p a t - — i = 0, 



respectively. Note that in the second-class constraint (|47|) 
we split the momentum conjugated to a according to re- 
lation (|20|) . and hence, the second-class constraint |47|) 
results in identity (f2"2l . Furthermore, second-class identi- 
ties (l46t and (I47|) will become auspicious at the quantum 
level since they enclose important operator identities. 



One can develop further the constraint T 2 (|45|) by ex- 
pressing the velocities in terms of the momenta, a = 
~{N/a) [p„/(( 7 - l)a 2 H 2 + 2)] and i = -SlN/{<y - 
l)a 3 H 2 by using the second gauge ([39|) . Thus, a direct 
calculation on the first term in (|4"5|) yields 



p-X 



N 



a[(7- l)a 2 iJ 2 + 2] 
Nfl 



(j-l)a 3 H 2 



Pt +p a a. 



(48) 



Finally, after a lenghty but straightforward computation 
the constraint T 2 can be cast as 
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F2 = Nlp 2 a 



(7 - l)a 3 H 2 



Pi 



Pad 

N 



a H 



J_ A72TT2 _ J_ 

or or 



n 



We could also use Eq. (|18p and the second-class con- 
dition (|47p . which impose the identities pt = — O and 
p a = 2aa/N, respectively, to conclude that the constraint 
T 2 ~ gives rise to a quadratic expression for the in- 
volved momenta. That is, the second gauge condition 
P9P shifts the problem from the linear dependence in 
the momenta, to a convenient quadratic expression for 
the physical momenta. To close this section we must 
mention that the constraints T\ and J-2 form an algebra 
isomorphic to the Lie algebra associated to the lower tri- 
angular subgroup of SX(2, K) as argued in the Appendix. 
This is the starting point to achieve an algebraic quanti- 
zation as we will sketch below. 



VI. QUANTIZATION 

In this section we study the canonical quantization of 
our system. Also, we will sketch an alternative differ- 
ent quantum theory for our model which emerge from 
the corresponding first-class symmetries. To this end, we 
emphasize the totally dissimilar nature which first- and 
second-class constraints play in the quantum theory, and 
also, we explore the different senses in which the physical 
states of our theory can be defined. 

We start in the conventional way by promoting the 
classical constraints into operators, densely defined on a 
common domain in a proper Hilbert space. As it is well 
known, we can only achieve a consistent classical theory 

I 



T 2 



we will work on the assumption that the commutators 
of these quantum constraints form a closed Lie algebra 
which will be also isomorphic to the algebra g. In fact, 
the classical first-class constraints are isomorphic to the 
algebra g associated to the lower triangular subgroup G of 
5L(2,R) (see the Appendix). Quantization of the lower 
triangular subgroup of SL(2, R) by algebraic methods 
was extensively studied in [27| (see also [28[ for com- 
parison). Now we explore the rather different senses in 



[( 7 -l)a 2 if 2 + 2]|. (49) 

I 

by implementation of the Dirac bracket. Once this is 
done, the second-class constraints are eliminated off the 
theory by converting them into strong identities. At the 
quantum level this is mirrored by defining the quantum 
commutator of two quantum operators as 

[A,B] :=i{XB}\ (50) 

where the Dirac bracket {•,•}* is defined in Eq. (|A7|1 . 
Thus, with this prescription the operators corresponding 
to second-class constraints are also enforced as operator 
identities [2q |. For our system, this yields the quantum 
operator expressions 

Si = fl v -^±=0, (51) 

S 2 =p a -^ = 0, (52) 

which, in particular, tell us the character of the quantum 
operators Il v and p a . For the rest of the variables, we 
choose to work on the "position" representation, where 
we consider the position operators by multiplication and 
their associated momenta operators by — i times the cor- 
responding derivative operator when applied on states 
defined on a suitable Hilbert space. 

By defining the quantum first-class constraints as 



(53) 



(54) 

I 

which the quantum constraints can be used to define ap- 
propriate physical states. 



A. Naive Dirac constraints 

First, we explore the Wheeler-DeWitt equation emerg- 
ing by considering the physical states ^ of the theory as 



:= -iN 




d 



(7 - l)a 3 H 2 dt 
xa [(7 - l)H 2 a 2 + 2] } , 
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FIG. 1: WDW potential for geodetic gravity with naive Dirac 
constraints. 




for the Dirac conditions on physical states. Let {C a } 
be a set of quantum constraints operators that gener- 
ate a nonunimodular gauge group with the commutators 
[C a ,Cb] — ifabCc, where /£, are the structure constants 
of the corresponding Lie algebra. Thus, the "unimodu- 
larization" procedure for nonunimodular groups dictates 
the consideration of the physical states \^f) as those satis- 
fying C a \^) = -(«/ 2 )/abl*)- Such modified Dirac condi- 
tions agree with the naive Dirac constraints if, and only 
if, the group is unimodular. 

Accordingly, for our theory the modified Dirac condi- 
tions for the gauge group invariant quantization of the 
system can be shown to be equivalent to 



those defined by naive Dirac conditions 

A* = o, 

= 0. 



(55) 
(56) 



Equation ([53)1 simply tells us that our physical states ^> 
are not explicitly depending on the phase space variable 
N. However, due mainly to the complexity of our WDW 
Eq. (|56p . we have not succeed in finding explicit solu- 
tions for the physically admisible quantum states. We 
note that the last term in the operator ([51} will bring a 
vanishing contribution to the WDW equation, and also 
we see that the ^-dependence can be avoided by assuming 



*(o,t) := 
equation 



-iQt 



if), where ip := ip(a) satisfies the WDW 



d 2 



#0 = o, 



where the potential U(a) is given by 



U(a) 



[(7-l)tfV + 2] 2 (l- 7 ffV 



(57) 



(58) 



which is recognized as the potential function found in [9] 
by repeated use of the master evolution constraint ([26} 
in Eq. ([S"6"} . The behavior of this potential is drawn in 
Fig. [l] where we can see the characteristic potential bar- 
rier. As discussed by Davidson and coworkers, it can 
be shown that after considering appropriate boundary 
conditions the big-bang singularity in our quantum the- 
ory can be neutralized by properly choosing the origin 
as inaccessible to wave packets. For further details on 
the behavior of the potential U (a) , the reader is referred 
to 1. 



B. Modified Dirac constraints 

As discussed in Refs. [HI, [30[, there exists a proce- 
dure which allows us to reduce nonunimodular groups 
to unimodular ones and this in turn brings a remark- 
able alteration for systems amenable to geometric or 
algebraic quantization which comprises a modification 



2J I*) 



0, 
0, 



(59) 
(60) 



which consequently define physical states |'F). Equa- 
tion (|59[) is equivalent to the homogeneity condition 
\\ft(rN)) = r-^^iN)) for r > 27]. Further, ([59} 
can be explicity solved by taking \^) = \ip), where 
A is a constant, and is a function of the variables a 
and t. Once more we also do not have control on the 
explicit solutions for the physically admisible quantum 
states. The ^-dependence can be avoided by assuming 
\ip(a,t)) := e~' int \tp), where \ip) is thought of as a func- 
tion of the scale factor a only, which satisfies the WDW 
equation 



da 2 



Ma) 



( 7 - l)a 2 H 2 + 2 



4a 2 



\<p(a))=0, (61) 



where the potential U(a) was described in the previous 
subsection. Hence, we see that our modified quantum 
theory brings out an extra potential term into our WDW 
equation, which succinctly differs from the one found 
with the naive Dirac procedure. 

We note that the extra term is purely emerging from 
the modified quantum Dirac Eqs. ([B"9"]) and l[6"0"]) , and it 
is completely absent while considering the naive Dirac 
procedure. This term will be nonvanishing even in the 
Einstein limit (7 — > 1), where it goes as aT 2 . Further 
studies about the possible physical implications of this 
term could be carried out. The behavior of the modi- 
fied potential is drawn in Fig. [21 where we can notably 
see that the central barrier potential present in Fig. [1] 
is almost vanishing while an infinite barrier emerges at 
the origin. Until our knowledge, the resulting unbounded 
potential is not realistic despite the first-class constraints 
suggest this modified description. Nevertheless, one can 
not resist the speculation of such possible quantum be- 
havior. Thus, rather than a nice potential, this time 
it is a more complicated function with distinct features 
notwithstanding the internal constraint symmetries that 
demand an unimodularization procedure. 
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FIG. 2: WDW potential for geodetic gravity with modified 
Dirac constraints. 



VII. CONCLUDING REMARKS 

By making use of the Ostrogradski formalism we have 
developed an alternative Hamiltonian description of the 
RT brane model. Unlike the Hamiltonian treatment by 
Davidson and coworkers for this model [13, [H| > our anal- 
ysis above keeps the original variables without the ne- 
cessity of introducing nondynamical variables. At first 
sight, this may look like an unnecessary complication 
since the configuration space is initially increased only 
to be reduced again at a later stage by imposing the con- 
straints and fixing the gauge. Nevertheless, it is hoped 
that despite computational complications we have pro- 
vided an improvement of physical clarity, in particular, 
concerning the geometrical meaning of the constraints 
of the theory and the physical content of the achieved 
canonical transformation. Although the Ostrogradski ap- 
proach has a price to pay, since neither the momentum 
P has the meaning of mechanical momentum nor Hq has 
to do with the energy of the system, as it is custom- 
ary, these quantities are adequate for providing a set of 
canonical equations which correctly describe the evolu- 
tion of the system. Also, an important point to mention 
is that the formalism is rich enough to demonstrate the 
real role of the extra terms coming from the surface: the 
phase space constraints of the system impose identities 
for these quantities which are valid at both classical and 
quantum levels, hence eliminating the unphysical degrees 
of freedom. 

In spite of the fact that this model is a second-order 
derivative theory and since it is well known that at 
quantum level the energy for higher-order derivative La- 
grangians is almost always unbounded below, the RT 
model results an exception due to its linear dependence 
on the accelerations which in turn contain important 
physical information commonly absent in higher-order 
derivative theories [3l|. Like it or not, until now a Hamil- 
tonian approach for the RT field theory demands the 
use of extra unphysical degrees of freedom at the be- 
ginning which by means of the phase space constraints 
are frozen out. Our adopted treatment renders the pas- 
sage to a full quantization for the system which can be 



achieved by means of an inspired canonical transforma- 
tion. We conclude further that the Ostrogradski quan- 
tum approach has exactly the same unique degree of free- 
dom as the Davidson and coworkers approach. Although 
our Wheeler-DeWitt equation for the scale factor is not 
analytically manageable, it is good enough to substract 
from it some interesting features. In particular, the po- 
tential we found is exactly the same as the one exten- 
sively discussed by Davidson et al. [9j. Furthermore, our 
Hamiltonian approach makes feasible the quantum treat- 
ment of Lagrangians with linear higher-order derivative 
dependence in the fields. 

It is suggesting that relativistic theories linear in the 
accelerations, for which characteristic surface terms are 
commonly neglected, are, as a generic feature, reluctant 
to quantization. To present day, quantization for these 
kind of systems have been mainly studied by considering 
some extra degrees of freedom by several other methods. 
From this point of view, our intention has also been to 
introduce the Hamiltonian-Ostrogradski approach as a 
geometrical powerful method to beset this sort of sys- 
tems. An specific example of this would be to apply 
our treatment to the almost forgotten idea concerning 
the rigid bubble electron, for which a linear correction in 
the extrinsic curvature of the electron surface is added to 
the Dirac-Nambu-Goto action in contrast to the conven- 
tional first order method where a surface term is omit- 
ted [331 . It will be worked elsewhere. 
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APPENDIX A: ALGEBRAIC PROPERTIES OF 
THE CONSTRAINTS 

We can construct the matrix Cab whose elements are 
the Poisson brackets of all the constraints Ca where 
A,B= 1,2,3,4. Hence, 



(Cab) 





f° 








1 











aN 













\0 


a 2 $ 


-N 3 n 




(Al) 



in the constraint surface. This matrix has rank 2, which 
is a signal that we have two first-class constraints [26| . It 
is also important to mention that constraints C\ and C3 
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form an algebra, namely, 



by 



{C 1 ,C 1 } = 0, 

{Ci,C 3 } = — c 3 . 
{C 3 ,C 3 } = 0, 



(A2) 



which reflects the invariance under reparametrizations of 
the RT field theory as a fundamental gauge symmetry. 
Indeed, this algebra results an isomorphism of the Lie 
algebra g associated to the lower triangular subgroup of 
SX(2, R) with positive diagonal elements, G. Such Lie 
algebra g is spanned by the matrices [33j 



h := 



whose commutator is 











U -J 


, e := 


! 


! 



[h, 



-2e~ 



(A3) 



(A4) 



and hence we realize the isomorphism through the iden- 
tification C\ i— ► h/2 and C3 1— ► e~. Among the relevant 
properties of the subgroup G we refer that G is two- 
dimensional, non-Abelian, connected, and nonunimodu- 
lar. This last property plays an important role in our 
quantum theory, as developed in section IVH 

Also, among the second-class constraints, (|3"r?)) and 
([37|) . we can construct the matrix Su = {Si, Sj}, given 



(Su) 



a$ / -1 



N 



and its inverse 



(0 = ^- 



N 
a<I> 



1 



1 

-1 



(A5) 



(A6) 



where I, J = 1,2. The matrix S IJ will help us to con- 
struct a Dirac bracket in the standard way: Let / and g 
be two arbitrary functions then 

{/, 9V := {/, 9} - Si}S ij {Sj, 5} , (A7) 

where { • , • } stands for the Poisson bracket defined in (|3ip . 
As it is well known, classically, the Dirac bracket is es- 
sential to eliminate the second-class constraint off the 
theory by converting them into simple functional identi- 
ties. The need for the Dirac bracket is also very relevant 
at the quantum theory since we can only reach a consis- 
tent quantization procedure through the implementation 
of this bracket $2$. 
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